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METHOD AND SYSTEM FOR MEASURING OPTICAL 
SCATTERING CHARACTERISTICS 

CROSS REFERENCE TO RELATED APPLICATION 

[0001] This patent application is a divisional of, and claims priority to, commonly 

owned United States patent application serial number 09/840,060, filed April 23, 2001, 
entitled, "Method and System for Measuring Optical Scattering Charateristics," by Leif 
Fredin, Robert Chin and William Hallidy, which will issue as United States Patent 
Number 6,606,148 on August 12, 2003, and which is incorporated herein by reference for 
all purposes. 

TECHNICAL FIELD OF THE DISCLOSURE 

[0002] The present disclosure relates in general to optical system monitoring, and, 

more particularly, to a method and system for measuring optical scattering 
characteristics. 

BACKGROUND 

[0003] Optical fibers increasingly constitute the chief means for transmitting 

information through the world's telecommunications network. Certain characteristics of 
an optical fiber can also be used to generate information rather than just transmit it. 
Specifically, the temperature of an optical fiber affects the amount and wavelength of 
light that will be scattered in response to a transmitted pulse. Careful measurements of 
scattered light can therefore be used to determine the temperature at points along an 
optical fiber. As another example, mechanical stresses on the fiber affect the amount of 
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certain wavelengths of light that will be scattered in response to a transmitted pulse. 
Once again, measurements of scattered light can yield useful information. 
[0004] Other optical systems also scatter light in correlation with characteristics 

of interest. For example, an air-filled region may scatter light in proportion to the density 
of pollutants or another constituent element of interest. Accurately measuring the extent 
to which certain wavelengths or ranges of wavelengths of light are scattered provides 
information about other characteristics of the system. 

[0005] In a conventional method a time-limited pulse of light with an 

electromagnetic spectrum of average wavelength X is produced at an excitation source 
and sent through an optical fiber. When the excitation source is a laser, the 
electromagnetic spectrum is often very narrow and is referred to in shorthand as a single 
wavelength. As the pulse traverses the fiber, backward scattered light is produced. Three 
types of backward scattered light, among others, are of interest: Stokes light, anti-Stokes 
light, and Rayleigh light. Stokes and anti-Stokes light are collectively referred to as 
Raman light. Stokes light constitutes an electromagnetic spectrum having an average 
wavelength greater than X. Anti-Stokes light constitutes an electromagnetic spectrum 
having an average wavelength less than X. Rayleigh light has the same wavelength X as 
the excitation source. The width of the Stokes and anti-Stokes spectra, as measured by 
the difference in wavelength between the points of 50% intensity, is often much greater 
than the width of the time-limited pulse spectra and the Rayleigh light spectra, especially 
if that pulse is produced by a laser. 
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[0006] Some of the Rayleigh, Stokes, and anti-Stokes light travels to the end of 

the fiber at which the pulse was introduced, while some is scattered at an angle such that 
it is absorbed by the cladding of the fiber or escapes. The location from which the 
backward scattered light originated can be determined by the time between the 
introduction of the pulse and the receipt of the light. After a pulse is introduced into the 
fiber, backward scattered light is continuously received and time functions of the total 
intensity across the Stokes and anti-Stokes spectra can be determined. Under particular 
circumstances, the temperature of a point in the fiber has a known relationship to the ratio 
of the anti-Stokes light produced at that point to the Stokes light produced at that point. 
If, however, the intensity of the excitation per area of the fiber core is too high, non-linear 
distortions eliminate the temperature proportionality. Increasing the measurement 
accuracy of Stokes and anti-Stokes intensity as a function of time without introducing 
non-linear distortion, increases the accuracy of the resulting calculation of temperature as 
a function of position in the fiber. 

[0007] United States Patent No. 5,113,277 discloses a Fiber Optic Distributed 

Temperature Sensor System. The '277 patent contemplates introducing a light pulse from 
a light source into a fiber. The scattered light is then divided by wavelength spectra with 
detectors positioned to receive the Stokes light and anti-Stokes light, respectively. The 
measurements made by the detectors are then introduced into an equation to determine 
the temperature at each measured distance. 

[0008] The use of timed pulses of light to detect temperature or mechanical stress 

can require expensive components. For example, a light source that has sufficient power 
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and produces light of a wavelength that has scattering characteristics allowing for 
measurements of scattering over a long distance of fiber can be very expensive. 
Additionally, the electronics necessary to convert the received intensity of back scattered 
radiation into a digital representation become more expensive as their processing speed 
increases. Increasing the spatial resolution of the temperature measurements using timed 
light pulses requires digital representations of back scattered radiation intensity for 
smaller periods of time. Such representations are only available with the use of faster, 
and consequently, more expensive electronics. Additionally, high power pulses can 
cause stimulated emission of Raman light. Such stimulated emission cannot be 
distinguished from backscattered radiation and renders calculations inaccurate. 
[0009] The time pulse method disclosed in the '277 patent also uses optical 

components to screen Rayleigh scattered light from the sensors. Analyzing the 
characteristics of Rayleigh scattered light can result in useful information indicating 
possible mechanical stresses in the optical fiber. This information is not available when 
the wavelengths comprising the Rayleigh scattering are blocked from the sensors. 

SUMMARY OF THE INVENTION 

[0010] A method and system of measuring optical scattering characteristics is 

disclosed. None of the advantages disclosed, by itself, is critical or necessary to the 
disclosure. 

[001 1] A system is disclosed for measuring optical scattering characteristics that 

includes a laser that produces an excitation signal. An optical fiber is coupled to the 
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laser. At least a portion of the excitation signal enters the optical fiber as a coupled 
excitation signal with a continuous waveform and an amplitude modulated at variable 
frequencies. A first detector receives radiation backscattered from the coupled excitation 
signal by the optical fiber. In a more specific embodiment, the coupled excitation signal 
has a power less than 500mW. In another more specific embodiment, the optical fiber is 
a single mode optical fiber. 

[0012] A method is disclosed for measuring optical fiber characteristics that 

includes coupling a continuous wave laser excitation signal to an optical fiber. Radiation 
backscattered by the optical fiber in response to the coupled excitation signal is detected 
to produce a backscattered radiation signal. The backscattered radiation signal is mixed 
with the excitation signal to produce a mixed signal. The mixed signal is filtered to 
reduce the magnitude of frequencies other than conjugate mixing frequencies relative to 
the conjugate mixing frequencies. The filtered signal is digitized and the magnitude of 
backscattered radiation from a specific portion of the fiber is calculated based on the 
digitized signal. In a more specific embodiment, the temperature of a specific portion of 
the fiber is determined from the magnitude of the backscattered radiation. 
[0013] It is a technical advantage of the disclosed methods and systems that 

backscattered radiation from an optical target receiving a variable frequency interrogation 
signal is detected. 

[0014] It is also a technical advantage of the disclosed methods and systems that 

less expensive electronics can be used to monitor variable frequency backscattering. 
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[0015] Another technical advantage of the system and method disclosed is that 

lower cost excitation sources producing less power can be used to produce accurate 
results. 

[0016] Another technical advantage of the system and method disclosed is that a 

lower power excitation signal can be coupled to a fiber to reduce non-linear distortion. 
[0017] Another technical advantage of the system and method disclosed is that 

the temperature at a specific point of the optical target can be determined. 
[0018] Another technical advantage of the system and method disclosed is that 

mechanical stresses of an optical fiber can be determined. 

[0019] Another technical advantage of the system and method disclosed is that 

the particle density of an air- filled region can be determined. 

[0020] Another technical advantage of the system and method disclosed is that 

the frequency difference between the detected reference signal and the detected 
backscattered radiation can be used to determine the origin of the backscattered radiation. 
[0021] Another technical advantage of the system and method disclosed is that 

the backscattering characteristics of single mode fiber can be determined. 
[0022] Other technical advantages of the present disclosure will be readily 

apparent to one skilled in the art from the following figures, descriptions, and claims. 
Various embodiments of the invention obtain only a subset of the advantages set forth. 
No one advantage is critical to the invention. For example, one embodiment of the 
present invention may only provide the advantage of detecting backscattered radiation, 
while other embodiments may provide several of the advantages. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0023] A more complete understanding of the present disclosure and advantages 

thereof may be acquired by referring to the following description taken in conjunction 
with the accompanying drawings, in which like reference numbers indicate like features, 
and wherein: 

[0024] FIGURE 1 is a graph of electromagnetic spectra; 

[0025] FIGURE 2 is a graph of a chirped, variable frequency, modulating signal; 

[0026] FIGURE 3 is a diagram of a system for measuring optical scattering 

characteristics in accordance with one embodiment of the present invention; 
[0027] FIGURE 4 is a graph of amplitude modulation frequencies of 

backscattered radiation; 

[0028] FIGURE 5 is a graph of a stepped, variable frequency, modulating signal; 

and 

[0029] FIGURE 6 is a diagram of a system for measuring characteristics of 

optical fiber in accordance with one embodiment of the present invention. 
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DETAILED DESCRIPTION OF THE DISCLOSURE 

[0030] FIGURE 1 is a graph of electromagnetic spectra. One point in the fiber 

will reflect some light of the same optical frequency as the light being transmitted, known 
as Rayleigh scattering 10. Of less intensity are reflections of light at optical wavelengths 
both longer and shorter than the incident light. The shorter wavelength light 12 is known 
as anti-Stokes scattering. The longer wavelength light 14, is known as Stokes scattering. 
The Stokes and anti-Stokes light are collectively known as Raman light. The intensity of 
each type of scattered light is a function of the intensity of the incident light. If the 
incident light is amplitude modulated, the amplitude of the backscattered radiation will 
also be affected. The different types of backscattered radiation are also affected by 
different characteristics of the fiber. For example, the amount of Rayleigh light scattered 
at a point in the fiber is related to the mechanical stress of the fiber at that point. The 
ratio of anti-Stokes light to Stokes light scattered at a point in the fiber is related to the 
temperature of the fiber at that point. 

[0031] Figure 2 is a graph of a chirped, variable frequency, modulating signal. 

The graph shows the frequency of the signal as a function of time. The modulating signal 
20 progresses linearly from a minimum frequency 22 to a maximum frequency 24 during 
a chirp period 30. The slope of the chirp 26 can be determined by dividing the change in 
frequency (the difference between the maximum 24 and the minimum 22) by the chirp 
period 30 or duration. The duration of the chirp 30 is preferably longer than the amount 
of time it takes for light to make a round trip through the fiber or optical target to be 
interrogated. The greater the proportion of the chirp period 30 to the round-trip time, the 
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greater the proportion of the detected radiation that can be used to measure optical 
backscattering characteristics. 

[0032] Figure 3 is a diagram of a system for measuring optical scattering 

characteristics in accordance with one embodiment of the present invention. A first 
frequency generator 52 provides a frequency chirped signal in accordance with Figure 2. 
The signal is amplified by an amplifier 80 and controls the output amplitude of a laser 50. 
In alternative embodiments the output of the laser can be directed to an external 
modulator that is driven by the first signal The output of laser 50 is a laser excitation 
signal. A low power laser excitation signal can be used to decrease nonlinear 
backscattered radiation responses. For example, a power of less than 500mW allows for 
use of a 1541 nm laser source while reducing nonlinear response. Alternatively, a higher 
power laser excitation signal could be reduced to less than 500mW by an external 
modulator or other optical device. The amplitude modulated light is directed to an 
optical target 75. In one embodiment, an optical fiber 74 is used to direct the amplitude 
modulated light to the optical target 75. In one embodiment, an optical fiber is the optical 
target. The optical fiber can be multi-mode fiber or single-mode fiber. A portion of the 
backscattered light from the optical target traverses filter 90 and is received by a 
avalanche photodiode 58. The filter 90 determines the type of backscattered radiation 
received by the avalanche photodiode 58. For example, the filter 90 may allow only 
Rayleigh radiation or Stokes radiation depending on the wavelengths that the filter 90 
transfers and blocks. 
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[0033] The avalanche photodiode 58 outputs a signal corresponding to the energy 

of photons received. That signal is amplified by an amplifier 80. A mixer 84 receives 
the amplified signal together with the modulating signal from the first frequency 
generator 52. The mixer 84 is a device that produces output signals at the sum and 
difference frequencies of the input signals. The mixer output signal is provided to a low 
pass filter 60, an analog-to-digital converter (ADC) 62, and a fast fourier transform (FFT) 
circuit 64. In an alternate embodiment, the FFT circuit 64 could be replaced by a 
software-implemented fast fourier operation. The digitized frequency information is then 
received by a processor 86. 

[0034] Figure 4 is a graph of amplitude modulation frequencies of backscattered 

radiation. A range of backscattered radiation frequencies are received at any particular 
time. The frequencies range from the minimum chirp frequency 22 to the maximum 
chirp frequency 24. The received frequencies are periodic over the chirp period 30. For 
example, at a time ti a range of frequencies 23 are received. The highest frequency was 
reflected by the nearest point in the optical target and is essentially identical to the current 
frequency of laser modulation. The lowest frequency was reflected by the furthest point 
in the optical target and is equal to the laser modulation frequency at a previous time. 
The time difference in modulation is equal to the time required for the light to traverse 
the optical target twice. For example, if the optical target is an optical fiber of length L 
and light travels through the optical fiber at speed c, the lowest frequency of received 
radiation will be the frequency at which the laser was modulated at a time previous. 
[0035] An expression for the instantaneous frequency of the chirp is: 
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/(/)=/ 0 +ymod(/ 5 T) Eqn. 1 

where f 0 is the minimum frequency, y is the chirp rate, and □ is the chirp period. The 
laser's output power waveform then has the form: 

P(t) = 0, /<0 



P(t) = P 0 {l-msin 2 [<p(t)/2]}, 0<t EqU ' 2 

with 

q>(t)=2njf(t)dt Eqn. 3 

where (p(t) is the phase of the waveform. The excitation shown in Eqn. 2 can be rewritten 
as: 

P(t) = P 0 {l-±m + -^cos[cp(0]} . Eqn. 4 



The time span of the dashed regions in Figure 4 is just *f . In one embodiment, only data 
from outside the dashed regions is considered for determining optical backscattering 
characteristics. In that embodiment, data is available if < x . 

[0036] Eqn. 4 shows that the intensity modulation of the backscattered radiation 

received from the optical target will have both a DC and an AC component. The AC 
component of the modulation of the received backscattered radiation intensity as a 
function of time from a fiber of length, L, and absorption coefficient, a (£) , can be 

expressed as an integral over the length of the fiber, after an initial transient period of one 
round trip time on the fiber, 2L/c: 
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r(t) = ±mP 0 jexp -l\a(V)dt 



Eqn. 5 



for 



t > 



where a(/) measures the returned strength, from position /, of the backscattered signal 
that is trapped in the fiber and c is the speed of light in the fiber. In one embodiment, 
a (i) is assumed to be a constant, independent of /, so that the interior integral of Eqn. 5 
is equal to □/. With that assumption the complex return, R(t) can be defined as 



Eqn. 6 



so that: 



r(0 = Re[R(0]. 



Eqn. 7 



[0037] Figure 3 illustrates that the real return signal is mixed with the real 

excitation signal at mixer 84. Because these signals each consist of both an AC and a DC 
component and because the AC part of each real signal is half the sum of the 
corresponding complex signal with its conjugate, the real mixed signal contains four 
types of terms: DC terms from the DC/DC mixing, terms at the original chirp frequencies 
from the AC/DC mixing, terms at twice the original chirp frequencies from the direct AC 
mixing (eV and eV 1 ) and low frequency terms from the conjugate mixing (eV 1 and e'V). 
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[0038] In one embodiment, the DC terms are eliminated by coupling only AC 

from the frequency generator 52 and avalanche photodiode 58 to the mixer 84, for 
example using a capacitor. The low pass filter 60 eliminates the AC/DC mixing terms 
and the direct AC mixing terms. The only terms passed in this embodiment to the ADC 
62 are the low frequency terms that result from conjugate mixing. 
[0039] The result of the mixing and filtering may then be written as: 



M(0 = Ae'* W R(0 



Eqn. 8 



Thus 



M(0 = ^je" 2a V(^)exp 



-iln J f(V)dt % 

A 2t 



dl 



Eqn. 9 



where K and A are constants that depend on circuit parameters. 

2L 

[0040] In one embodiment we restrict the acquisition time to — < t <x , resulting 

c 



in: 



J /(*')*' = [/. +Y*]| 



2£ 
c 



2t\ 1 (It 



Eqn. 10 



Writing / (t ) = f 0 + y t and substituting k (t) = ^£^1 9 we have 

c 



M(t)mM (k) = A)e- 2ae *(t)e^ e nm dl . 



Eqn. 11 
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where the dependence of k on t is implicit. Performing a Fourier transform of M{k) 
yields an expression that can be solved for a (/) . 

00 

M (/)ee $M(k)e i2nk! dk 

Eqn. 12 

= Ae- 2ai a(i)e^~) , 0<1<L 

2L 

However, data is available only over the finite range of time, — <t<x, and over 

c 

corresponding ranges of f{t) and k(t), where 0<^ <k<k 2 . Thus, for a continuous 
range ofk, we have 

M\l,k x ,k 2 }^ jM(k)e i2nkl dk Eqn. 13 

*i 

which, after some manipulation, yields 



sinc[(* 2 -*,)(/ 



Eqn. 14 



where (k) = —(k l +k 2 ). If the remaining terms under the integral were slowly varying, 
the sine function could be approximated by a Dirac delta function. However, for positive 
k, we note that (k) > -^(k 2 - k x ) , so that this condition is not met in practice. 
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[0041] We are only able to collect data for positive frequencies, / (t) , so we only 

have experimental data for positive values of k(t). However, if we examine Eqn. 11 for 
M (k) , we see that if the condition 

□yDL/c) 2 «l Eqn. 15 

is met, then M(-k) «M*(fc), where * indicates complex conjugation. Thus, M{k) 
may be extended by this process, which we call bookmatching, to include negative k. 
This will permit us to set (&) = 0 in Eqn. 14. The remaining terms will then vary much 

more slowly than the sine function. If we also choose ^=0, Eqn. 14 takes the 
approximate form 

M~(i 9 k 2 )= \M(k)e i2Kki dk 

* 2Ak 2 \e- 2ae a(£)smc[2k 2 ('-')]<*' Eqn. 16 

o 

*Ae- 2a 'a(T), 0<1<L 

where we have approximated e ^ «land 2& 2 sinc [2^ (/-/)]«8 (/-/). The first 
approximation becomes exact when y = 0 . 

[0042] When backscattering from the fiber occurs from both discrete and 

continuous scatterers, the received signal from different locations on the fiber may 
include a wide dynamic range. In that instance, the contribution in Eqn. 16 from the 
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sidelobes of the sine function for a discrete scatterer may swamp the distributed 
scattering signal from nearby locations. In order to reduce this effect, we multiply the 
mixed signal of Eqn. 1 1 by a low pass window function, W(k) , to produce 

M w (k) = W(k)M(k) Eqn. 17 

By defining W(k) as a real, symmetric function of k, the function M w {k) will satisfy the 
condition for bookmatching if M{k) does. 

[0043] After bookmatching, we follow the procedures for M w {k) which led to 

Eqn. 16 for M{k) to obtain 

M ~w (' > K ) = } M w (k)e i2n "dk 

I 

*A$e- 2ae cj(£)w(i-ey£ Eqn. 18 

0 

*Ae- 2ai a{iy 0<1<L 

where W (l) is the Fourier transform of a suitably normalized window function, W(k), 
chosen to ensure that W ) approximates the Dirac function. 

[0044] Typical, non-normalized window functions that we have used include the 

offset Gauss type 

W G (k) = l, 0<k<k 0 

= e {kl ' ko) , k 0 <k<k 2 Eqn. 19 

- 0, elsewhere 
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and a generalized raised cosine function 
W c (k) = l, 



( 


71 


(k 






Y 


COS 












V 


2 


k 






) 



= 0, 



0<fc<£ 0 

, K<k<k 2 

elsewhere 



Eqn. 20 



where ko and a are parameters that determine the width and rate of decay of the window. 



mrM<n 



function, Wg, may be written 



m*)=g«-iim 

*2 



Eqn. 21 



where G(k) is a Gaussian function and Y\(k) is the unit step function with support 
-k 2 <k<k 2 . The Fourier transform of this widow is given by 

W G (e) = G(£)®2k 2 - sine (2k 2 £) Eqn. 22 

where G(^) is also a Gaussian function and ® denotes a convolution. 
[0046] The central peak of a sine function is twice as wide as each of its 

sidelobes. We can choose G(k) in Eqn. 21 so that the width of G(i) matches that of 
the sine's central peak. From Eqn. 22, this allows us to significantly reduce the sidelobes 
of W G (£) in comparison with those of the sine function while minimizing the spread of 
its central peak. 
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[0047] Data will normally be collected in discrete samples at equal intervals 

rather than continuously. Furthermore, it is also desirable to avoid taking zero frequency 
(DC) data. By choosing sample times so that the least sample frequency is half the 
frequency interval between samples, we arrive, after bookmatching, with a set of equally 



spaced samples, |m [k n ) : -N < n < N - 1 j , where k n = 



frequency interval 6/ =y8t and sample time interval 5t , 



n + — |5A:with Sk=^^- 
2 c 



[0048] 



Discrete Fourier transform of the data yields 



M(l,k N )=f d M(k n )e i2 « k " i Sk 



n=-N 



= ASk\e 2at a(e) 



sin[^(/-^)] 



dt 



Eqn. 23 



2f 

where k N = ^ . The function in curly brackets is called the array factor. It will 
c 

approximate a Dirac delta function over the region of integration provided the following 
conditions are met: 



2Uk N L»\ 



Eqn. 24 



and 



UUkL«\ 



Eqn. 25 



If these constraints are satisfied, we find 
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M {i,k N )*Ae- 2ai a(j), 0<1<L 



Eqn. 26 



in agreement with our earlier result from Eqn. 16, where the right hand side Eqn. 26 is 
independent of k N . If we relax the assumption that a (£) is a constant, independent of/, 
that was made in simplifying Eqn. 5 to Eqn. 6, we find 



M ,4 exp -2ja(£)d£ a(/), 0<1<L Eqn. 27 



[0049] Thus, we see that if the constraints of Eqn. 15, Eqn. 24 and Eqn. 25 

are satisfied, the transformed, low pass part of the mixed signal allows us to determine 

the signal backscattered from fiber as a function of position, / , along the fiber. We also 

observe the expected exponential decay of the return signal with / . 

[0050] Figure 5 is a graph of a stepped, variable frequency, modulating signal. 

Like the chirped signal 20, the frequency of the stepped signal 40 varies with time. 

Unlike the chirped signal 20, the stepped signal 40 varies discretely, not continuously. 

The stepped signal 40 progresses from a minimum frequency 42 through discrete 

frequency steps to a maximum frequency 44. Each step has a set duration 46. 

[0051] Figure 6 is a diagram of a system for measuring characteristics of optical 

fiber in accordance with one embodiment of the present invention. A first frequency 

generator 52 is shown producing a first signal that directly modulates a laser 50. In 

alternative embodiments the output of the laser can be directed to an external modulator 

that is driven by the first signal. The first signal has a frequency that varies with time. In 
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one embodiment the signal of Figure 3 is used to modulate the amplitude of the laser 50. 
A splitter 70 reflects a reference signal and transmits an interrogation signal. The 
interrogation signal is directed into an entrance 72 coupled to an optical fiber 74. 
Radiation backscattered from the fiber 74 exits 76 the fiber. In one embodiment a PIN 
diode 56 converts the reference signal to a detected reference signal. In an alternative 
embodiment, the detected reference signal is replaced by a copy of the first signal. The 
detected reference signal, like the first signal, is electronic rather than optical. An 
avalanche photo diode 58 converts the backscattered radiation into a detected 
backscattered signal. Each of the detected signals are coupled to an amplifier 80. 
[0052] A second frequency generator 54 is phase-locked to the first frequency 

generator 52 such that a second signal is produced that has a fixed frequency offset from 
the first signal. In one embodiment that fixed frequency offset is 10 kilohertz. The first 
mixer 84 receives the amplified detected backscattered signal as one input and the second 
signal as another. It outputs a first mixed signal. The second mixer 82 receives the 
amplified detected reference signal, or first signal copy, as one input and the second 
signal as another. It outputs a second mixed signal. The first and second mixed signals 
are coupled to a low pass filter 60 to remove high frequency components. The mixed 
signals are then digitized by analog-to-digital converters 62. The digitized signals are fed 
to fast fourier circuits 64 to be transformed to the frequency domain. In an alternate 
embodiment, the fast fourier operation could be accomplished in software rather than 
hardware. The offset frequency domain signal, 10 kilohertz in this embodiment, of each 
of the digitized frequency domain signals is captured in the peak capture circuit 66. The 
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offset frequency contains the modulation information. A divider circuit 68 divides the 
captured offset frequency of the first mixed signal by the captured offset frequency of the 
second mixed signal and provides the result to a processor 86. In an alternate 
embodiment, the captured offset frequency of the first mixed signal can be divided by the 
captured offset frequency of the second mixed signal in software. The processor 86 
determines characteristics of the optical fiber 74 from the received signals as expressed in 
the divider output. 

[0053] Taking y = 0 in the chirped frequency response of Eqn. 1 1 yields an 

expression for the complex mixed signal of the form 

L 

M(k) = A je^'cj (£) e i2nkt d£ . Eqn. 28 

0 

In this case, the conditions of Eqn. 15 for bookmatching are met exactly. If data is 
collected by stepping the laser modulation through equally spaced frequencies, and if the 
conditions of Eqn. 24 and Eqn. 25 are met, then the discrete Fourier transform of the 
sampled, mixed signal, yields Eqn. 26 or Eqn. 27. 

[0054] Although the present disclosure has been described in detail, it should be 

understood that various changes, substitutions and alterations can be made thereto 
without departing from the spirit and scope of the invention as defined by the appended 
claims. 
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